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ABSTRACT 

A technique Is presented for calculating feed- 
back and feedforward gain matrices that enable a VTOL 
aircraft to track input commands of forward and ver- 
tical velocity while maintaining acceptable responses 
to pilot inputs. Leverrier's algorithm is used for 
determining a set of state-variable, feedback gains 
that force the closed-loop poles and zeros of one 
pilot-input transfer function' to pre-selected posi- 
tions In the s-plane. This set of feedback gains is 
then used to calculate the feedback and feedforward 
gains for the velocity-command controller. The 
method Is computationally attractive since the gains 
are determined by solving systems of linear, simul- 
taneous equations. The method has been used In a 
digital simulation of the CH-47 helicopter to control 
longitudinal dynamics. 
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INTRODUCTION 

In recent years the overcrowding of this country's 
major airports has led to Increased Interest In the 
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development of commercial VTOL'or iSTOL' aircraft 
capable, of operating in the 0-500 mlle'irmge. These 
vehicles would, operate from separate runways at 
existing airports or from rooftops- or short runways 
in or near business districts. Studies have shown 
(1^ , (2) that such aircraft could reduce the total 
trip time providing they do not have to operate 
under existing take-off and landing procedures with 
their long delays. This will require innovations in 
vehicle design plus Improved navigation, guidance, 
and control systems. One approach to this latter 
problem Is to use on-board digital computers to 
handle the navigation, guidance, and control functions 
in an adaptive mode. Thus it is necessary to have a 
controller capable of following guidance commands from 
the computer and at the same time present acceptable 
flying qualities to pilot inputs. 

The procedure to be presented requires a set of 
state-variable feedback gains that result In accept- 
able responses to pilot inputs. This set of gains la 
then used to generate feedback and feedforward gains 
that maintain the same response to pilot Inputs, but 
enable the vehicle to track input commands from the 
guidance system in the form of changes in vertical 
and forward velocity. The procedure has been used 
in a digital simulation of the longitudinal dynamics 
of the CH-47 helicopter. 

BASIC CONCEPTS 

In considering the longitudinal dynamics of the 
CH-47 helicopter, the linearized equations of motion 
may be written In the form 


Ax + G6 


( 1 ) 


where 


u (forward velocity perturbation) 
w (vertical velocity perturbation) 
q (pitch rate) 

6 (pitch angle perturbation) 


6 (differential collective perturbation) 
e 

6^ (collective perturbation) 


■ 

) 


A is a 4 X 4 coefficient matrix 
G is a 4 X 2 control matrix 


Introducing state-variable feedback from all 
states^ to both inputs results in a control vector 
of the form 
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« - Kx + « (2) 

— — — p 

where 


K Is a 2 X 4 matrix of feedback gains 
5 Is a 2 vector of pilot commands 
-P 

Applying this control to equation (1) yields 
X - (A + GK)x + (3) 

The Laplace transform of equation (3) Is 


[si - (A + GK)]x(s) - G« (s) (4) 

~ ~P 

This equation leads to eight closed-loop transfer 
functions relating the four state variables to the two 
controls. These transfer functions may be obtained by 
employing Cramer's rule. Thus, they all have the same 
denominator, |sl - (A + GK) | , but different numerators 
that are obtained by taking the determinants of the 
matrices that result when the appropriate columns of 
the G matrix are substituted for columns of 
(si - (A + GK) ] . Each transfer function Is a ratio 
of polynomials In s and the coefficients In these 
polynomials are functions of K. If values of K can 
be found to force these coefficients to take on pre- 
determined values, then the poles and zeros of the 
closed-loop transfer functions can be placed any- 
where In the s-plane. 

DETERMINATION OF FEEDBACK GAINS 


While It la not possible to specify all of the 
coefficients of all transfer functions, several 
techniques (2)»(4),(^ have been presented for cal- 
culating the feedback gains necessary to place the 
poles of all transfer functions and the zeros of 
one pre-selected function at arbitrary locations In 
the s-plane. The methods referenced above are. all 
based on Leverrler's algorithm and the techniques 
of (4) and compute the gains by solving linear 

simultaneous equations. 

Before the gain matrices for the velocity con- 
troller can be calculated, feedback gains must be 
determined that yield acceptable responses to pilot 
Inputs. To accomplish this, the procedure developed 
In (50 was employed. This method Is based on the 
fact that the closed loop numerator and denominator 
functions of this system may be written as 


*l'^e. 


J"1 J c 
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( 6 ) 


2 4 4 4 X. ,6 

^ ^ k N(x ,6 )+ I I 1 ( 7 ) 

^ 1-1 J-1 J ^ 1-1 J-1 *j'^c 


where N(x^,6^) Is the open-loop, numerator of the 
transfer function and D-. is the open-loop denomina- 
tor. The coefficients of these terms may be generated 
using Leverrler's algorithm. 




represents the determinant of [si - A] with 


i column replaced by and the J column replaced 
by (columns of the control matrix G) . The coeffi- 
cients of these terms do not come directly from 
Leverrler's algorithm, but can be determined using 
the procedures of either (^ or (2) . 

Equations (S) and (6) show that the numerator 
polynomials are linear functions of (n - 1) of the 
2n feedback gains. Thus specifying the coefficients 
of one numerator allows direct calculation of (n - 1) 
elements of one row of the gain matrix. Arbitrarily 
picking the n^^ element of that row and substituting 
Into (7) allows direct calculation of the remaining 
n elements of K. 

One should note from the above discussion that. 

In general, the poles and zeros of only one transfer 
function may be arbitrarily placed and that the gain 
matrix Is not unique. 

SELECTION OF A MODEL 

While, in theory, the above technique Is capable 
of generating feedback gains that can place the 
closed-loop poles and zeros of one transfer function 
at arbitrary locations In the a-plane. In practice, 
considerable care must be exercised In specifying 
these values. If the pole-zero locations or the 
arbitrary gain element' are not carefully selected, 
either the gains will be too large or the other trans- 
fer functions of the system will have zero locations 
that result In unacceptable responses. For the CH-47 
the poles were placed to satisfy settling time require- 
ments on step Inputs. However, little or no informa- 
tion was available on how to specify zero locations 
and the method employed was basically trial and error. 

CLOSED-LOOP RESPONSE TO PILOT INPUTS 

Once the elements of the gain matrix have been 
determined so that the poles and zeros of one transfer 
function are fixed, the numerators of the remaining 
closed-loop transfer functions must be checked. While 
all transfer functions have the same denominator so 
that If one Is stable all will be stable, stability 
is not the only factor to be considered In evaluating 
step responses. Whether the steady-state change to 
a unit step Input Is too large or too small must be 
considered and right-half-plane zeros close to the 
origin may cause problems (7) . 

The closed loop numerator and denominator coeffi- 
cients are easily obtained using Leverrler's algorithm 
with A replaced by (A + GK) . With the coefficients 
known, any polynomial root-finding routine can be 
used to determine the pole and zero locations, and 
any transfer function with a right-half-plane zero 
should be checked by plotting its step response. 

If the step Input responses of all transfer 
functions are acceptable i then a satisfactory pilot- 
command control system of the form shown below has 
been generated. 



Fig. 1 Pilot-Command Controller 
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CLOSED-LOOP VELOCITY COMMAND CONTROLLER 

With the above system It Is possible, using the 
steady state responses of the transfer functions, to 
calculate the control Inputs required to produce 
command changes In forward and vertical velocity 
(u and w ). However, this would be an open-loop 
vefoclty controller with the Inherent disadvantages 
of open-loop systems. To overcome this, the following 
system configuration Is proposed. 



Each of the above equations represents a third order 
polynomial In s with four coefficients. While It Is 
not possible to control all of the coefficients. It 
la possible to control one coefficient from each equa- 
tion. If the system Is to track velocity commands. 


the steady-state response of u and w to a unit step 
Input on u must be one and zero respectively and the 
steady-staEe response of u and w to a unit step Input 
on w must be zero and one respectively. This requires 
that^^the constant term of the numerator of the u/u 
transfer function equal the constant term of the ^ 
denominator (which is known) and the constant term of 
the w/u transfer function equal zero to provide 
decoupling. Placing these constraints on equations 
(13) and (14) yields two linear simultaneous equations 
to solve for k_ and k_ . Similar restrictions on 
2ll 221 

the constant terms of the w transfer functions result 
in two equations that yleld'Ht, and k, ■. Once the 


‘12 


‘22 


values of are known can be determined using 
K - [K^ - K2’] (17) 


SIMULATION 


Fig. 2 Configuration for Velocity-Command Controller 
where 


10 0 0 
0 10 0 


Is a 2 X 2 matrix of feed-forward gains 

Is a 2 X 4 matrix of feedback gains 

c. Is a 2 vector of commanded changes In u and w 

jc* “ Is a 2 vector of the perturbations in u and w 

From figure 2 It can be seen that 


This design procedure has been applied to the 
linearized longitudinal dynamics of the CH-47 heli- 
copter. The results to be presented are for nominal 
velocities of ISO kts forward and 250 fpm vertical 
descent. The A and G matrices shown below are based 
on stability and control derivatives provided by 
NASA's Langley Research Center. 


A “ 


.05191 -.03898 8.8944 -32.176 

.02731 -.57793 250.67 1.2324 

.00014 .01769 -1.3044 0 

0 0 1.0 0 


± - K^x + K2 [c - K^x] 
j5_ - Kj^x - + K^c 

- " ^*^1 ■ *^ 2 '^- * h- 

(8) 

(9) 

(10) 

G - 

-.14909 

.01857 

.31973 

0 

-1.2698 
-8.9842 
.22782 
0 . 


or 


The 

open-loop 

characteristic equation for this A 

_6 - + K2£ 

(11) 

matrix is 
L 

1 

1 

Thus the closed-loop dynamics are defined by 

(12) 

®0L 

The 

- 8 + 1.934s- - 3.579s - .221s + .0117 

sign changes In the coefficients Indicate an 

X “ (A + GK)x + GK 2 £ 

unstable system and right-half-plane poles occur 


Now let K be the set of feedback gains obtained by 
specifying the poles and zeros of one of the pilot- 
command transfer functions. Then It Is apparent that 
the velocity-command and pilot-command transfer func- 
tions will have the same poles. Employing Cramer’s 
rule and expanding the numerator determinants It can 
be shown that the numerators of the u/u , w/u , u/w , 
and w/w transfer functions may be expressed is ** 


N(u,u ) 
c 


N(w,u ) 
c 


N(u.w^) 


N(w,w^) 










(13) 

(14) 

(15) 

(16) 


at .034 and 1.2. 

Specifying the poles and zeros of the w/6 
transfer function as 

Poles •> -.75, -.8, -.8 ± j.A 
Zeros “ -1.0, -.8 ± j.4 

resulted In the following gain matrix: 

-5. 


K - 


.0667 

.0021 


-.02 

. 0034 


-23.75 

28.08 


>.17 . 
.324_ 


when specified at -.02. This feedback gain 

matrix'*' yields denominator and numerator polynomials 
for pilot-command transfer functions as shown below 
with the steady-state response of each listed to the 
right. 
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Steady-state 
Response /Inch 


D..,- 

s%3.15s-’+3 

.88s +2.2S+.48 


of control 

CL 

N(5 ,u)- 

.15s^ - 9 

.238^ 

-16 

. Is -6 

.33 

-13.2 

fps 

e 

N(6 ,w)“ 

.02s^ - 

.6s^ 

-s 


.51 

- 1.06 

fps 

e 

N(6 ,q)- 

C 

.32s^ + 

.216^ 

+ 

.01s 


0 


N(6 ,0)- 


.328^ 

+ 

.21s + 

.01 

.02 

rad 


N(6^,u)- -1.27s^ -10. 4s^ 

N(6 ,w)- -8.98s^ -23. 4s^ 

® 3 2 

N(6 ,q)- .228s + .0l3s - 

2 

N(6 ,e)- .228s + 

c 


-11.4s -1.58 

-21.6s -7.19 

.006s 

.013s- .006 


fps 

fps 


- 3.29 
15.0 

0 

- .013 rad 


Thus , the system has been stabilized and the 
steady-state responses look reasonable. Sign changes 
In the u/6 , w/5 , and 8/6 numerators Indicate right- 
half-plane®zerosf however, in u/6 and w/6 they are 
located at 64 and 34 respectively and do not seriously 
affect the responses. The 6/6 right -half-plane zero 
is at .14 and does cause some ^ase reversal. However, 
6 is not very responsive to 6 commands and this may 
not be objectionable. 

To determine the feedforward matrix, equations 
(13) - (16) are used to form the following 2x2 
system of equations. 


-6.33 


- .51 


and 


-6.33 
- .51 


-1.58 


-7.19 


-1.58 


-7.19 





.48 


^11 

- 





0 


21 J 


. 


*1 

k, 


0 


2 i2 

« 



k- 


.48 


^22 




The resulting feedback and feedforward matrices 

are 


-.0105 

-.003 

-23.75 

-5.17 

.0034 

-.065 

28.08 

.324 

-.0772 

.017'] 



.0054 

-.068 




These gains, used in Che system configuration of 
figure (2) , yield the following veloclty-comnand 
transfer function numerators and steady-state responses: 


Steady-state Re- 
sponse to Unit 
Step Command 


N(u,u )- 

-.018s^+ .66s^ +1.17S +.48 

1.0 


c 

N(w,u^)- 

-.050s^- .08s^ - .04s +0 

0.0 


N(q,u^)“ 

-.023s^- .02s^ - .0008s 

0.0 


N(e.u^)- 

- .023s^- .02s -.0008 

.002 

rad 

N(u,w )- 

.089s^+ .553a^+ .503a + 0 

0.0 


c 

N(w,w )- 

.611s^+1.58a^ +1.45S +.48 

1.0 


N(q.w^)- 

-.Ola^ + .003s^+ .0006a 

0.0 


H(0,w )- 
c 

- .Ola^ + .003a +.0006 

0 

0 

rad 


T 

The responses to velocity commands of c^ - [10,0] 
and c_^ ~ [0,10] were plotted and ace shown in figures 
(3) and (4) respectively. 



U 

fps 



Tune, Sec 




[ 0 , 10 ] 


The Velocities settle to within 5Z of the desired 
steady-state changes in 5 seconds with no overshoot 
and the decoupling between forward and vertical velo- 
city commands is good. The pitch angle to vertical 
velocity command response shows a phase reversal that 
reaches a peak of 0.9 degrees in 1.3 seconds, but 
this may not be objectionable. If these responses 
were unacceptable to either pilots or passengers, 

Che gains could easily be recalculated using new 
pole-zero specifications. 


CONCLOSIONS 


A procedure has been presented chat allows direct 
calculation of feedback and feedforward matrices for a 
closed-loop velocity-command controller without affect- ' 
Ing the pilot-command transfer functions. The method 
has a very attractive computational characteristic in 
that it does not require a non-linear or iterative 
search technique to calculate Che gains. All of the 
gains are determined by solving linear simultaneous 
equations. Numerical results are presented for Che 
linearized longitudinal dynamics of the CH-47 helicopter 
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